Stability of multi-field cosmological solutions in the presence of a fluid 
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We explore the stability properties of multi-field solutions in the presence of a perfect fluid, as 
appropriate to assisted quintessence scenarios. We show that the stability condition for multiple 
fields <f>i in identical potentials Vi is simply d 2 Vi/dcj>i > 0, exactly as in the absence of a fluid. A 
possible new instability associated with the fluid is shown not to arise in situations of cosmological 
interest. 
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I. INTRODUCTION 

In models with multiple scalar fields, those fields may 
act collectively to drive an accelerated expansion even 
if each field individually were unable to, a phenomenon 
known as assisted inflation This possibility has now 
been widely explored for early Universe inflation, and 
has begun to be considered for quintessence models of 
the present acceleration as well 0t3]- ^ n * ne simplest 
scenario, a number of fields with identical uncoupled po- 
tentials may be invoked. 

Such scenarios clearly allow for solutions where the 
fields evolve together, by symmetry. This is not sufficient 
however to demonstrate assisted behaviour, as one must 
also check that such solutions are stable. In cosmologies 
where the fields are the only components, this was proven 
for the original exponential potential case in Refs. [2, [![, 
and Ref. [9( subsequently provided the general condition 
for stability which is simply that the potential for the 
fields be convex, d 2 Vi/dcf> 2 > where cf>i are the fields 
and each Vi(4>i) has the same functional form. 

In this article, we extend the result of Ref. @ to allow 
for the presence of a fluid, hence obtaining the stability 
conditions appropriate for assisted quintessence scenar- 
ios. 



of state of the fluid, and dots represent derivatives with 
respect to synchronous time t. 

A powerful method of exploring the stability of such a 
model is to move into a Hamilton- Jacobi type formalism 
[Tol [TT| in which the role of the clock is invested in a mat- 
ter field so that one can get rid of the unphysical degree 
of freedom represented by shifts in the time coordinate t. 
Here we choose one of the scalar fields <fio = cj> as the new 
time coordinate and perturb only the other n — 1 fields 
4>i. The actual field chosen is arbitrary; the instability 
we are exploring is in the difference between fields and 
would be identified whichever is chosen as the reference 
field. Primes will indicate derivatives with respect to the 
field. The Klein-Gordon equation for <f> is then a con- 
straint determining the relation between different time 
coordinates. To rewrite the equations of motion (JXJ) , @ 
and §3§ in this new formalism, first we combine them to 
obtain 
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where 



n-l 



(4) 



(5) 



II. BACKGROUND EVOLUTION 



p' — p/4>, H' — H j 4> and here and in what follows i runs 
from 1 to n — 1. Rewriting Eq. ([1} as 



n-l 



We consider a flat Friedmann-Robertson- Walker back- 
ground with n scalar fields and a fluid. The equations of 
motion are then the Friedmann equation, fluid equation, 
and n Klein-Gordon equations 



n-l 



3H 2 = p + W + -Y i - 



• 2 

h ; 



i=0 

p + 3Hjp = ; 
3H<p + diW = 0. 
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where reduced Planck units Mp\ = c = 1 have been used, 
H is the Hubble parameter, W = W((j)o, 4> n -i) is 
the potential with both self-interaction and interaction 
terms, 7 = 1+p/pisa constant giving the equation 



3iJ 2 = p + W- 



(0) 



1=1 



then substituting in Eq. Q , we get our Hamilton- Jacobi 
equation 



2H' - 3AH 2 + A(W + p)=Q. (7) 
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Note that if we send p, p' — > 0, we recover the Hamilton- 
Jacobi equation used in Ref. [9(. 

In the next section we will be writing down an equation 
of the form SX' = MSX where X = (H,p,(f> i ,(f>' i ) T is a 
2n-vector, so we combine Eq. Q with Eq. ^ to obtain 



9jAH 2 p - QHH'p' = 0. 



(8) 
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Finally we use the fact that 

d 2 = dSd 4 ) = -(3ff<£ + W')d$ + <j> 2 d 2 , (9) 

to get n — 1 equations of the form 

$4>'l - W'tf* + diW = 0. (10) 

Note that this equation as written is independent of H 
and p, though in practice those quantities will influence 
the evolution of 4>. 

III. STABILITY 

As in Ref. [|[ we wish to discuss the classical stability 
of solutions of the form 



<f>i(t) = (j){t) i = 0,...,n-l, 



(11) 



such that the fields evolve together. Solutions of this 
form are only possible if 



d,W\ 



= v\4>) v*,j, 



(12) 



where primes represent derivatives with respect to <f> and 
V is defined by this equation. If the fields are mutually 
decoupled, then each can be written with its own poten- 
tial Vi obeying 



Vi(x) + Ai = V(x) V %. 



(13) 



Here the A, are constants and can all be absorbed into 
some A = acting as a cosmological constant, so 

that W = J2i ^ + A. We will consider A to be negli- 
gibly small or zero. The co-evolving field solutions can 
also exist within a limited class of models with cross- 
couplings between fields; for further discussion on this 
class of solutions and the types of potential that permit 
it, see Ref. @. 

We are now in a position to apply linear perturbation 
theory to analyse the stability around solutions of the 
form Eq. ([IT]). When perturbing Eqs. Q, © and ([ID)) 
around Eq. dill) , some simplifications arise since 4>[ = 
and A — n on the background (8 A — 2 £\ <5<^)- Perturb- 
ing the Hamilton-Jacobi equation we find 



SH ' = ASH + BSp + CS& + D<5(/>- , 



where 
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From perturbing Eq. ([8]) we obtain 

Sp' = E6H + FSp + G64>i + H£$, 



where 
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QCHp' 
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(Note that we use Roman typeface to label coefficients 
and that in particular H, the Hubble parameter, and H, 
the coefficient of dfy, should be distinguished.) 
Finally, from perturbing Eq. (|10[) we obtain 



li54>i + JSfti, 



where 
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We thus have a 2n x 2n matrix of the form 
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such that, suppressing the subscript is where possible, 
the system is described by 



/6H'\ 
Sp' 



/A B C D 

E F G H 

1 

,0 L J 




(27) 



To analyse the solutions we use the characteristic equa- 
tion \M — A/ 2n |=0 to find eigenvalues. We then require 
that all eigenvalues satisfy 3?(A) < for stability. By 
standard manipulation of the determinant, the charac- 
teristic equation can be reduced to the block-diagonal 
form 



A- A B 

E F — A 
-A 1 
I, J- A 



= 0, 



(28) 
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and so the terms C, D, G, and H do not influence sta- 
bility. This expression demonstrates that perturbations 
in the scalar field are independent of perturbations in 
the metric 5H and the fluid Sp. In this paper, we shall 
concentrate on the case where the fields are decoupled 
such that djdiW{<t>) = V" djW = and j = i 
in Eq. ([25]) . but since the scalar fields' perturbations are 
unaffected by the presence of a fluid, the more complex 
models discussed in Ref. Q are still valid here. 
The eigenvalues are 



A 4 



and 



H.sp = \ (A + F ± v/A 2 - 2AF + 4EB + F 2 ) (29) 



J ± y/P 



41 



(30) 



Eq. (|30p is a (n — 2)-degenerate pair. We use the con- 
vention <j) > and V < as we are interested in situa- 
tions where the field is rolling down the potential. Con- 
sequently, if we substitute Eq. (I2"5j) back into Eq. (13TJ1) . 
taking account of the simplification above 



V 
20 2 



1± 



A<j) 2 V" 



y/2 



(31) 



we find that the requirement for stability is simply 

V" > . (32) 

This is intuitively the case as we know from being able 
to write the determinant in block-diagonal form, that 
the scalar fields evolve independently of the metric and 
thus evolve independently of each other. So for the per- 
turbations to die away we require that the fields further 
down the potential are evolving more slowly than those 
further up the potential, namely that the potential is 
convex. This is unchanged from the result of Ref. @, i.e. 
the stability within the scalar sector is unchanged by the 
presence of the fluid. 

Analysing Eq. (I2TJ1) is less straightforward. In the ab- 
sence of a fluid the requirement for stability is simply 
that found in Ref. M, 



T <0 ' 



(33) 



which is automatically satisfied. This can be derived by 
noting that in this limit the product BE vanishes, and F 
is irrelevant as there is no Sp, and so the RHS of Eq. (|2"9")l 
simply becomes A which can be simplified to the above 
form. In the presence of the fluid, there is a new pertur- 
bation degree of freedom corresponding to a shift in the 
fluid density with respect to the value of the scalar field 
4> acting as the clock, which could in principle yield a 
new instability. As stability requires that the real part of 
both eigenvalues is negative, we need to check that both 
A + F < and BE - AF < are satisfied. 



The condition for the former to be satisfied can be 
shown, after some algebra, to be 

n<j) 2 36ff 3 n( 7 + I)0 2 + 6H~/ 2 p (p - 6H 2 ) + jp4> 

^ 7 2 p(6ff 2 + p) (&H 1P + ^ , (34) 

provided 

7P - n<j) 2 ^ 0. (35) 

For general 7 the condition for the latter is rather long, 
so we show only the simpler case of matter domination 
(7 = 1): 



l08H 4 n 2 ^ 3 § p (QHp + , 
provided 



ruf>[3H + <f>] +p 



p-nc/) 2 ^ 0. 



(36) 



(37) 



Both conditions have a direction of inequality that de- 
pends on a subsidiary condition. We first note that these 
will not be satisfied in complete generality, because when 
the subsidiary condition is saturated the main condition 
is not, meaning there must be some regime where the 
main condition fails in the vicinity of the region where 
the subsidiary condition flips. 

However, in situations of practical cosmological inter- 
est, Eqs. (|35|) and (|37|) will both be satisfied with a pos- 
itive inequality, indicating the fluid dominating over the 
combined field kinetic energy. This is of course true in 
the early stages of cosmological evolution where the fluid 
dominates the total density, but it remains true up to 
the present, even though the total field energy density 
dominates, as the observed equation of state w = 7 — I 
close to —1 requires that the field kinetic energy is still 
subdominant to the fluid. Only sometime in the future 
is the field kinetic energy expected to overtake the fluid. 

Qualitatively, taking the kinetic part of the scalar field 
to be subdominant, <f> 2 <C p + W, in Eq. (f3~4"|) we require 
LHS < RHS which is satisfied thanks to the LHS be- 
ing suppressed by a factor of (f) 2 . Similarly the LHS of 
Eq. (|36|) is suppressed by a factor (ft 3 . 1 

We can make this argument quantitative by defining a 
parameter a that measures the dominance of the fluid: 



p = 3H 2 (l + a) 



n<t> 2 



W 



6H 2 3H 2 



(38) 



1 We also note that were we in a kinetic-dominated case, revers- 
ing the required inequality, those suppressions would become en- 
hancements suggesting that stability would be restored deep in 
the kinetic-dominated regime. 
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Reintroducing Planck masses for clarity, then Eq. ([36 
becomes 

108#V</> 3 < 162H 5 M^(a + l) 3 
+162H 4 M^n(a + l) 2 </> 



+MH 3 Mpw(a + 1 



,212 



9H 2 M^(a + l) 2 



+9HM$ l n(a + l)<f> 2 + 3M^n(a + 1) 



(39) 



where each term is positive. 

The first situation of interest is when we have matter 
domination with the scalar fields highly subdominant, as 
it is interesting to know if we have already reached the 
solution (ITT1) by the time the Universe starts to accel- 
erate. Since a — > in this situation and we have that 
n<f) 2 <C 6H 2 Mp it we see that even the first term on the 
RHS is sufficient to guarantee the satisfaction of the in- 
equality. 

The second case we are interested in is our present 
one, in the vicinity of the onset of acceleration. Then the 
energy density of the scalar field and fluid are of similar 
order. The condition for acceleration is 



a > 



ne < 1 - 



IP 
2H 2 ' 



where 



\3H 



2H' 



2nH 2 



(40) 



(41) 



is a slow-roll parameter. Note that for a single scalar 
field and no fluid, we recover the first Hamilton-Jacobi 
slow-roll parameter 



(42) 



We thus see that stability in this situation is a little more 
complicated. If p > H 2 then stability is guaranteed but 
if p < H 2 then acceleration can take place for either 
sign in Eq. (|35|) and since equality does not occur at 
the same time in any, let alone all the conditions (|3"4")l 
(|35|) (|36]) . stability is no longer guaranteed. Having said 
that, we know at the present time the equation of state 
parameter w ~ — 1 and so we still expect Eq. (|3"5j) to be 
positive and thus, since 1 + a is still of order 1, Eq. (f39|) 
is again satisfied. Using the same approach one can show 
that Eq. (f3~4"|) is also satisfied under these conditions and 
hence we have stability both at early times and during 
conditions similar to the present day. 

We end this section with a brief comment on the inclu- 
sion of a separate constant A term. For the scalar field 
sector the situation is unchanged by the introduction of a 
fluid; Ref. Q showed that introduction of A could change 
the type of stability (oscillatory or non-oscillatory about 
the background solution) , but could not change stability 
to instability. However the new fluid degree of freedom 



obeys a more complex equation and it does not appear 
possible to make a general statement on the effect of A. 
Nevertheless, cosmologies with a fluid and fields are al- 
ready capable of explaining observations without needing 
A and in this set-up a zero value of A is the most natural 
choice. 



IV. CONCLUSIONS 

We have derived a set of stability conditions against 
classical perturbations for multi-field cosmological solu- 
tions in the presence of a fluid. Our focus has been on 
the co-evolving ('assisted') class of solutions for which a 
general stability analysis has been possible; for more gen- 
eral cross-coupled potentials such solutions typically will 
not exist and any stability properties must be assessed on 
a case-by-case basis. Within the scalar sector, the con- 
dition for stability is independent of the presence of the 
fluid, as seen immediately from the block-diagonal form 
of the characteristic equation, Eq. (|28l) . indicating that 
field perturbations are decoupled from the fluid sector. 

The introduction of the fluid introduces a new degree 
of freedom, corresponding to a shift in the fluid density 
with respect to the scalar field 'clock'. Stability under 
this perturbation does not appear to be guaranteed for 
arbitrary backgrounds, but we have shown that for situa- 
tions of physical interest this stability condition is always 
met. 

Our stability criteria are local, rather than global, con- 
ditions on the solutions; since they depend on the condi- 
tions on the potential they may be satisfied during some 
parts of the evolution and not others. An example would 
be multiple fields on a cosine potential as in axion mod- 
els, e.g. as in Ref. 12], where the co-evolving solution 
would be unstable around the maximum (V" < 0) and 
stable around the minimum (V" > 0). In regimes where 
the co-evolving solution is unstable it would clearly be 
extremely fine-tuned, whereas in those where it is stable 
there will be a basin of attraction around the co-evolving 
solution making co-evolving evolution much more likely. 

Our results also introduce a new perspective on track- 
ing solutions [ill, [lU . Normally one considers the fluid 
evolution p(a) as the fixed quantity; the field <p is then 
shown to have various possible behaviours as a function 
of a, which converge together if the tracking criteria are 
met. In our formalism, instead it is 4>i(4>) which is fixed 
to be equal to <j>, and p is perturbed as a function of 
<f>. It would be interesting to make a more extensive ex- 
ploration of tracking properties in the Hamilton-Jacobi 
formulation. 
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